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Abstract 

The relaxation times over which dissipative fluxes restore their steady state values have been 
evaluated for a pion gas using the 14-moment method. The effect of the medium has been im¬ 
plemented through a temperature dependent mr cross-section in the collision integral which is 
obtained by including one-loop self-energies in the propagators of the exchanged p and a mesons. 
To account for chemical freeze out in heavy ion collisions, a temperature dependent pion chemical 
potential has been introduced in the distribution function. The temperature dependence of the 
relaxation times for shear and bulk viscous flows as well as the heat flow is significantly affected. 
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I. INTRODUCTION 


Characterizing the thermodynamic properties of matter composed of strongly interacting 
particles has been the premier objective of heavy ion collision experiments at the Rela¬ 
tivistic Heavy Ion collider (RHIC) at Brookhaven and the Large Hadron Collider (LHC) at 
CERN l|. Relativistic hydrodynamics has played a very important role in analyzing the 
data from these collisions j2| and providing a viable description of the collective dynamics of 
the produced matter. Recently, the observation of a large elliptic flow(n 2 ) of hadrons in 200 
AGeV An-An collisions at RHIC could be explained quantitatively using a small but hnite 
value of shear viscosity over entropy density {rf/s) |3[. However, a consistent formulation 
of relativistic dissipative fluid dynamics is far from trivial. The hrst order theories are seen 
to lead to instabilities due to acausal propagation of perturbations. The second-order the¬ 
ory due to Israel and Stewart ^ currently appears to be the most consistent macroscopic 
formulation to study collective phenomena in heavy ion collisions BQ. 

Though the hydrodynamic equations may be derived from entropy considerations using 
the second law of thermodynamics, a microscopic approach is necessary in order to determine 
the parameters e.g. the coefficients of shear and bulk viscosity, thermal conductivity and the 
relaxation times of the corresponding fluxes. The Boltzmann transport equation has been 
used extensively as the underlying microscopic theory to estimate the transport coefficients 
of relativistic imperfect fluids. In this approach the (differential) scattering cross-section in 
the collision integral is the dynamical input and plays a significant role in determining the 
magnitude of the transport coefficients. The case of a pion gas has received some attention 
and several estimates of the transport coefficients exist in the literature. In all the cases the 
TTTT cross-section corresponds to the one in vacuum. Either the chiral Lagrangian has been 


used 

data 


7] to derive the scattering amplitude or it has been parametrized from phase shift 


12l |. In general, medium effects affect the collision integral in two competing ways. 


The larger phase space occupancy due to the Bose factors (1 -|- /)(! -|- /) for the hnal state 
pions results in an increase of the collision rate. This is compensated to some extent by a 


smaller effective cross-section on account of many-body effects (see e.g^ 
effect of the medium on the viscosities |l^ and thermal conductivity 15| was studied using 


13|). Recently, the 


the Chapman-Enskog approach and signihcant modihcation in the temperature dependence 
of the coefficients was observed. 
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In addition to the coefficients of viscosity and thermal conductivity the corresponding 


relaxation times r also go as input in the viscous hydrodynamic equations 


16 1. They 


indicate the time taken by the fluxes to relax to their steady state values and consequently 
play an important role in determining the space-time evolution of relativistic heavy ion 
collisions. This is more so for systems where r is of the same order or larger than the mean 
collision time tc of the particles since several collisions may occur during the relaxation of 
the dissipative flows to their steady state values as in the case of a strongly interacting 
system like the one created in heavy ion collisions. Moreover, though the magnitude of the 
shear viscosity is usually much larger than the bulk viscosity, the corresponding relaxation 
times may be comparable. Also, the ratios of the viscous coefficients to their relaxation 


times are found 


17| to behave differently with respect to temperature. There are a few 


estimates of the relaxation times available in t 
of the relaxation times have been evaluated in 


le literature. The temperature dependence 


10 


18 l | with a parametrized cross section 


which is independent of temperature. Constant values of transport coefficients have been 
used in and in jo] these quantities have been obtained using conformal held theory. 

In the present study we investigate the effect of the medium on the relaxation times of 
the dissipative hows. As is well known, the Chapman-Enskog approach leads to a linear 
relationship between the thermodynamic forces and the corresponding irreversible hows. 
Because of the parabolic nature of the equations of motion this results in inhnite speeds 
of these hows. In order to have access to the relaxation times we use the more general 
14-moment method due to Grad 19|. With the inclusion of the viscous pressure tensor 
and the heat how to the original (hve) hydrodynamic variables the relations between huxes 
and forces contain time derivatives of the huxes and cross-couplings between them. The 
hyperbolic nature of the equations of motion in this case result in hnite relaxation times of 
the dissipative hows. Our aim in this work is to estimate the change in the temperature 
dependence of the relaxation times for the shear and bulk viscous hows and the heat how 
for a hot pion gas on account of the in-medium cross-section. We thus evaluate the tttt 
scattering amplitude with an ehective Lagrangian in a thermal held theoretic framework 
and use it in the Uehling-Uhlenbeck collision integral which contains the Bose enhancement 
factors for the hnal state pions. In addition to a signihcant medium dependence we hnd 
the relaxation times for the viscous and heat hows for a chemically frozen pion gas to be of 
comparable magnitude. 
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The formalism to obtain the relaxation times of the dissipative fluxes using the 14-moment 
method is described in the next section. This is followed by a discussion on the medium 
dependent tttt cross-section in Sec.III. The results are given in Sec. IV and a summary in 
Sec.V. Details of calculations are given in Appendices A, B and C. 


II. THE RELAXATION TIMES IN THE 14-MOMENT METHOD 


We begin with the relativistic transport equation for the phase space density 


P^d^f{x,p) = C[f] 


( 1 ) 


where the collision integral is given by 


in which dTo = , 

^ {2'KfE 


C[f] = J dTkdTp, dTk'[f{x,p')f{x,k'){l + f{x,p)}{l + f{x,k)} 

-fix,p)f{x, k){l + f{x,p')}{l + f{x, k')}] W , (2) 

d?q 


, E = \/^ + and the term 


M-' = ~(2!r)V(p+«:-p'-A^') 

contains the differential scattering amplitude ■^. The factor 1/2 is to account for the 

dkl 

indistinguishability of the initial state particles which are pions in our case. In the moment 
method one attempts to obtain an approximate solution of the transport equation ([1]) by 
expanding the distribution function f{x,p) in momentum space around its local equilibrium 
value when the deviation from it is small. We write 


f{x,p) = f^^\x,p) + 5f{x,p), Sf{x,p) = f^^\x,p)[l + f^°\x,p)](j){x,p) 


( 3 ) 


where (j) is the deviation function. The local equilibrium distribution function is given by 


/(0)(a;,p) = 


exp{ 


p^Up,{x) - pi{x) 
T{x) 




-1 


( 4 ) 


in which T(a;), and p(x) are identified with the temperature, fluid four-velocity and 

pion chemical potential respectively. The latter arises on account of conservation of the 
number of pions due to chemical freeze-out in heavy ion collisions and is not associated with 
a conserved charge (see later). 
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Putting ([3]), the left hand side of ([T]) splits into a term containing derivative over the 
equilibrium distribution and another containing derivative over 0, 

= -£[ 0 ] ( 5 ) 

where the collision term reduces to 

^[0] = j dVk dVp, dTk'f^^\x, k){l + f^°\x,p')}{l + f^^\x, k')} 

[(j){x,p) + (j){x, k) — (j){x,p') — (j){x, k')] W . (6) 

To simplify the hrst term on the left hand side of eq. dSD the partial derivative over 
is decomposed with respect to the fluid four-velocity into a temporal and spatial part by 
writing = u^D + V^. Here D = is the convective time derivative and 
is the spatial gradient. The projection operator is dehned as where the 

metric = diag(l, —1, —1, —1). Using this, the local rest frame, where = (1,0), is 
dehned following Eckart as = 0. In this frame the spatial components of the particle 

four-current = J drpp^^f{p) vanishes. 

When the derivative over is taken with the above prescription we obtain a set of terms 
containing space and time derivatives over the thermodynamic quantities. While the space 
gradients of thermodynamic variables lead to thermodynamic forces, the time derivatives 
are eliminated using the equations of motion containing the dissipative huxes. These are 
discussed in Appendix-A. After some simplihcation we get 




X 


(r - - (n^n,)(v^n") + QV% - tU^Du^ 


T 


Tn 


+ r[{r(l - Y) + (7" - l)h - 7 '"} o - 


5 

P 


riT 




(7) 


with = p^/T, T = p ■ u/T and Q = Q/T'^, where, Q = -|- {p ■ — 7 '} + P ■ 

m{(7" — l)h — 7'"T}. The notation indicates a 

spacelike, symmetric and traceless form of the tensor The reduced enthalpy per particle 
is dehned as, h = h/T and P, n and stand for the pressure, particle density and heat 
how vector respectively. The 7 ’s and the (5’s are dehned in Appendix-A. 

For the remaining two terms in ([5]) we need to dehne the deviation function 0 and its 
derivative. Since the distribution function is a scalar depending on the particle momentum 
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and the space-time coordinate x^, the deviation fnnction is expressed as a snm of scalar 


prodncts of tensors formed from p^ and tensor fnnctions of x^. Following 


20| we write 0 as 


4){x,p) = A{x,t) - Bf,{x,T){Il^) + Cf,„{x,T){Il^Il'') 


( 8 ) 


where ( 11 ^) = and {fl^II^') = -|- A^"A^^ — lA^^A"^]!!^!!/?. 

Now the X and r-dependent coefficient fnnctions A, and are farther expanded in 
a power series in r snch that the last power is the one which gives a non-zero contribntion 
to the collision term, getting 

2 

A{x, t) = Ao + Ai{x)t + A 2 {x)t‘^ = ^ As{x)t% (9) 

s=0 

1 

B^{x, t) = Bof,{x) + Bir{x)T = '^{Bs)f,{x)T% (10) 

s=0 

C^u{x,t) = Co^u{x) ( 11 ) 

where in the last eqnation s = 0. This leaves ns with six x-dependent coefficients Aq, Ai, 
A 2 , Bqh, Bi^ and It is convenient to express them in terms of the thermodynamic 

flnxes (irreversible flows). 

Let us start with the viscous pressure which is dehned as, 

n = 11 £ir„A„„pV/<")(i + /(»>),#.. ( 12 ) 

Putting 0 from eq. ([ 8 ]) we get 

n = -T^A 2 J Tp (13) 


in which terms containing Aq and Ai vanish due to the properties of summation invariance. 
Note that only the scalar coefficients of (j) a ppe ar owing to the fact that only inner product 


of irreducible tensors of same rank survive 
we hnally get 


n = ?7.TQ!2^2 • 

We next turn to the energy 4-flow which is dehned as 


20|. Dehning J dTpf^^^l + /(°))gr- 

(14) 


P = / dTpp^A^^ip ■ u - h)/(°)(l + /(°))0. 


(15) 
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On putting 0 it retains only the vector coefficients in it by virtue of inner product properties 
of irreducible tensors and yields 

I, ^ j TpWA^^ir - h)T{W)f^^\l + 

= , (16) 

with/3„ = — f dTpf^^\l +—h)Afj,^p^^p'^. Finally we have the traceless viscous 
tensor dehned by 

{un = jvp (i7) 

Proceeding as before we dehne the tensor coefficient in 0 as 

(nn = -lp7o(C?‘'). (18) 

with 7 „ = —^ / dTpf^^\l + (ppPu){p^p'^) and p = m.„n is the mass density. 

To obtain the remaining coefficients Ao{x), Ai{x) and Biy{x) of equation ([9]) and (1T0|1 
we utilize the conservation laws obtained by asserting that the number density, energy 
density and the hydrodynamic 4-velocity can be completely determined by the equilibrium 
distribution function. This leads to the following constraint equations for the deviation 
function 0 


J dP,p%/W(l + /(°))(/. = 0, 

(19) 

dPp(p%)V(°)(l + /(°))0 = O, 

(20) 


(21) 


Putting the value of 0 in the eqs. (IT^ and (l2n|l we obtain relations involving coefficients 
A giving, 

< 1 iAq -|- (X2A1 -\- 03^.2 = 0, (22) 

CL 2 AQ -\- 03^1 -|- (X4^A2 = 0, (23) 

with an = f dPp/^°^(l -|- Since A 2 is known from (1TT|1 the other two A’s can be 

determined. Similarly the relation between the coefficients Bp coming from eq. dsn is, 

S°A^^ 6 o + BlA^^^bl = 0, (24) 
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with = J dFpf^^^l + /(°))r"'(n^)(n^). Details of the calculation are discussed in 

Appendix-B. 

Using equations fll4lll6lll8ll22ll23ll24p we obtain the complete set of coefficient functions in 


terms of the thermodynamic ffows. These are given by 

(0204 — al) n 


Aq = 
Ai = 

A 2 = 
Bqu = 
Biu = 

(cn = 


° (0103 —a|)nTa2 

(ciiCl 4 — 02^3) If 


0403) nTa2 


n 


nTa2 

nm ^ 60 ' 


A;<f 

<! 


nTPi 

5 

Pl« 




(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


Dehning all the space-time dependent coefficients appearing in eq. (IHD in terms of the known 
functions it is now possible to specify the deviation function 0 completely. We now use it 
in eq. ([5]) to evaluate the equations of motion for the dissipative fluxes. 

1. Bulk viscous pressure equation 


Taking inner product of both sides of eq. ([5]) with and applying the (inner product) 
properties of irreducible tensors 2^ we obtain the equation of motion for bulk viscous 
pressure, 


n = civ„«'‘ - ;^{/?c-Dn + acV^/pi, 


(31) 


where 


^ 1 T , an - 2a2a3ai + aial , 

Pc= —- 2 -+ 05 }, 

a^n 02 — oiOs 


(32) 


CXD 

The terms S'" appearing above are defined as S^{z) = E (x) denoting 

k=l 

the modified Bessel function of order I with z = m^/T. The integrals an, bn, otn and /3„ 


8 












have been defined earlier. Their explicit forms as well as those of the 7 ’s and 5’s expressed 
in terms of 5“ appear in Appendices D and A respectively. Also the number density n can 
be expressed as n = -^m^TS\. 

Retaining only the first term on the right hand side of (]3T|l the equation for the bulk 
viscous pressure reduces to the same as in the first order theory of dissipative fluids. The 
coefficient of this term corresponds to the bulk viscous coefficient C, and is given by 

C = T^. (34) 

0-22 

The quantity 022 in the denominator of fl34D stems from the collision integral containing the 
interaction cross section and is described in Appendix C. 

Note that eq. flMll contains the time derivative of the bulk viscous pressure and is hyper¬ 
bolic. This yields a relaxation time for bulk viscous flow given by 

= C Aft ■ (35) 


2. Heat flow equation 


In this case we take the inner product of both sides of equation ([5]) with {n^)r. Following 
similar techniques as above we get the equation for heat flow, 

1 


= TX 




T 


nh 


nT 


+ 7AV„(n'"') + aW'-n} 


(36) 


with 


/5a = 

7a = 


ax — 


Lili + 

/?i 7 o n 

3T 11 0204 — Uq , 

— IT —1^1-2 + 

n Pi 0/2 Oitts — 02 


C 1104 


“ 0203 , , 621 

-+ & 3 I +y]- 

oiOs h 


(37) 

(38) 

(39) 


The factors pi and /9i expressed in terms of Sf are given in Appendix- 

enthalpy per particle h can be expressed as h = 2 ;%. 

*^2 

The thermal conductivity in this approach is given by 

611 ’ 


D. The reduced 


( 40 ) 
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where hn follows from the collision term and is dehned in Appendix-C. 

So from the above equation the relaxation time for heat flow is obtained as 


ta = AT—/3a. 

nl 


( 41 ) 


3. Shear viscous pressure equation 

Multiplying both sides of equation (|5]) with (11^11'^) and using similar techniques as before 
produces the equation of motion for shear viscous pressure and is given by, 


(n^^) = 


2{W) - —{l3,D{Un - 


(42) 


with 


I3v = 




23 12 


6 rr / SI si,. 


(43) 


(44) 


The shear viscosity rj is dehned as the coefficient of the hrst term of (H2|l . In this approach 
it is obtained as 

Here again cqo contains the dynamics of interaction and is described in Appendix-C. 
From (l42ll the relaxation time for shear viscous how is given by 




(46) 


III. THE IN-MEDIUM Trvr CROSS-SECTION 


It is clear from the above expressions that the microscopic dynamics concerning the tttt in¬ 
teraction which governs the transport coefficients enters through the cross-section appearing 
in the collision integral. Adopting a phenomenological approach we consider tttt interaction 
to occur via a and p meson exchange using the well-known interaction 


^ = 9p^ ■ ^ X 


(47) 
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TT 


: vrvr scattering with self-energy corrections 


FIG. 1 


where n and a denote the isovector pion and scalar sigma helds respectively and the conplings 
are given by Qp = 6.05 and = 2.5. As indicated in fig. [1], the meson propagators in the 
s-channel diagrams are replaced with effective ones obtained by a Dyson-Schwinger snm of 
one loop self-energy diagrams involving the pion. The matrix-elements for vrvr scattering in 
the isospin basis is then given in terms of the Mandelstam variables s, t and u as 

s — t 


Mi=o — 2g 


s — u 

t — rn?p 


+ 


2 2 

J<7 <7 


u — m; 


+ 


+ 


-M/=i = gl 


s — + So- t — ml u — m 

2{t — u) ^ t — s 


s — + S 


u — mi 


u — s 
t — m^ 


2 2 


1 


1 


u — mi 


(48) 


t — ml 

where we have ignored the non-resonant 1 = 2 contribntion. The terms So- and Sp appearing 
in the isoscalar and isovector amplitndes respectively denote the vacnnm self-energies of the 
a and p involving only the pion in the loop diagrams. The cross-section obtained from 
the isospin averaged amplitnde |ATp = X]/l-^-fP/+ 1) is shown by the dashed 
line in fig. |2] and agrees very well with the estimate based on measnred phase-shifts given 
in {q]. In this way it is ensnred that the dynamical model nsed above is normalized against 
experimental data in vacnnm. 

In order to obtain the in-medinm cross-section the self-energy diagrams are now eval- 
nated at hnite temperatnre nsing the techniqnes of thermal held theory in the real-time 
formalism 


22 


23| . For the a meson only the ttyt loop diagram is calcnlated in the medinm 


whereas in case of the p meson in addition to the tttt loop graph, ttcj, tt /ii , vrai self-energy 
graphs are evalnated using interactions from chiral perturbation theory 2^. In contrast to 
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So- which is scalar, Sp has longitudinal and transverse parts and are dehned as 251 


s^ = -i(s;i + ^Soo), s^ = is 


2 ^ 00 ? ^00 — tie’ll 


(49) 


where is the 4-velocity of the thermal bath and (f = {u ■ q)^ — q^. The momentum 


dependence being weak 


25| we take an average over the polarizations, 

E, = i| 2 E^ + 4 r"E''] . 


(60) 


The imagina ry p art of the self-energy obtained by evaluating the loop diagrams can be 
expressed as [22| 

{ 2 Try Au^Uh 

- f^^\i^h))S{qo -u^-Uh) 

+ - f^°\uJh))S{qo -UJn + Uh)] + 

L2{{,f^^\uJh) - f‘'^\(^n))S{qo +Un- Uh) 

_(1 _ /(°)(a;^) - f^^\uh))S{qo + + Uh)}] (51) 

where is the BE distribution function with arguments cUtt = \Jk‘^ + rnl 


and Uh = y {q — kY + rnY- The terms Li and L 2 relate to factors coming from the vertex etc. 


details of which can be found in 


22l |. The angular integration is done using the ^-functions 


which dehne the kinematic domains for occurrence of scattering and decay processes leading 
to loss or gain of p (or a) mesons in the medium. In order to account for the substantial 
Stt and pn branching ratios of the unstable particles in the loop the self-energy function is 
convoluted with their spectral functions 26|, 

1 /■("ih+ 2 rh )2 


S(g,m/,) = 


N> 


dM^ X 


h J{mu-2VhY 


Tm 


TT -ml + iMVh{M) 


Y.(q,M) 


(52) 


with 


'‘{mu+2Th)'^ 


Nu = 


dM^ X 


2 (m^- 2 r ^)2 

—Im 


vr -ml + iMVh{M) 


(53) 
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FIG. 2: The vrvr cross-section as a function of Ec.m.- The dashed line corresponds to scattering in 
vacuum. Dot-dashed line refers to the in-medium cross-section involving only the pion loop for the 
a and p mesons. The solid line corresponds to the additional loops in the p meson self-energy. 


in which rrih and Th stand for the vacuum mass and width of the unstable meson h. The 
contribution from the loops with these unstable particles can thus be looked upon as multi- 
pion effects in tttt scattering. 

In relativistic heavy ion collisions, below the crossover temperature inelastic reactions 
cease and this leads to chemical freeze-out of hadrons. Since only elastic collisions occur the 
number-density gets fixed at this temperature and to conserve it a phenomenological chemi¬ 
cal potential is introduced which increases with decreasing temperature until kinetic freeze- 
out is reached . In this work we use the numerical results of the temperature-dependent 
pion chemical potential from the work of 28|] where the above scenario is implemented. It 
is depicted by the parametric form 


/i(T) = a + hT + dT^ 


(54) 


with a = 0.824, h = 3.04, c = —0.028, d = 6.05 x 10 ® and T, /i in MeV. 

We plot in £g. [2]the total tttt cross-section defined by (j(s) = - f dVL-— with —- = . 

2 d\l d\l 64:71^ s 

The increase in the imaginary part of the self-energy due to scattering and decay processes 
in the medium results in enlarged widths of the exchanged p and a. This is manifested in 
a suppression of the magnitude of the cross-section and a small downward shift in the peak 
position as a function of the c.m. energy. The effect is larger when more loops are included 
in the p self-energy compared to the pion loop as shown by the solid and dash-dotted lines 
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FIG. 3: Relaxation time of bulk viscous flow as a function of T. 


in fig. El 


IV. RESULTS AND DISCUSSION 


We now present the results of numerical evaluation of the relaxation times. We start with 
the relaxation time of the bulk viscous flow r^, as a function of temperature. In £g. [3] two 
sets of curves are displayed which correspond to two different values of the pion chemical 
potentials. The curves in the lower set (with circles) are evaluated with a constant value 
of pion chemical potential. We take this value to be /i = 85 MeV which is representative 
of the kinetic freeze-out in heavy ion collisions. The upper set consisting of curves with 
triangles show results for the temperature dependent pion chemical potential fi = fi{T). 
The three different curves in each set show the effect of the medium on account of the tttt 
cross section. The dotted curves show results where the vacuum cross-section is used. We 


have checked that our estimates for this case agree with 9|, [l^ for constant values of the 


pion chemical potential. The dashed curves depict medium effects corresponding to the pion 
loop in the a and p propagators. The relaxation times appear enhanced with respect to the 
vacuum ones indicating the effect of the thermal medium on Finally the uppermost solid 
curves correspond to the situation when the heavy mesons are included in the p propagator, 
i.e. for vrh loops where h = vr,a;,hi,ai. The larger effect of the medium on seen in 
this case is brought about by a larger suppression of the cross-section which appears in the 
denominator. The clear separation between the curves in each set displays a significant 
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0.1 0.12 0.14 0.16 
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FIG. 4: Relaxation time of heat flow as a function of T 


effect brought about by the medium dependence of the tttt cross section. Also the two sets 
of curves appear nicely separated showing the significant difference caused by a constant 
value of the pion chemical potential /i = 85 MeV and a temperature dependent one which 


decreases from 85 IV 


eV at kinetic freeze-out (T = 100 MeV) to zero at chemical freeze-out 


(T = 170 MeV) [291]. Consequently the upper set is seen to merge with the lower one at 
T = 100 MeV. 


Next we plot in fig. H] the relaxation time for the irreversible heat flow, t\ against tem¬ 
perature for the same two different values of pion chemical potentials mentioned above. In 
each set the curves are plotted for different ttyt cross sections. Similar to the earlier case 
here also we notice that the medium modified cross sections evaluated at hnite temperature 
influence the temperature dependence of t\ which appear enhanced for the in medium cases 
with respect to the vacuum ones. The multi-pion loop contribution due to heavier mesons 
in the p propagator turns out to be more significant than the ttyt loop. 


Finally we present the results for i.e, the relaxation time of the shear viscous flow in 
fig- 0 The curves are seen to follow the same trend corresponding to the cases described 


above though the magnitudes are a little lower. For t 


lower set of curves, this has already been seen 


in jg. 


le case of constant p, given by the 


10|. Recalling that Tr^ = ril3r^/nT, 


its smaller variation with temperature is due to the fact that the increase of p with T is 
largely compensated by the decrease in 1/nT, remaining approximately constant in the 
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• ■ • vacuum, |i=85 MeV 
•- • Tz-Ti, (1=85 MeV 
7i-h, (1=85 MeV 
A ■ * vacuum, |i=(i(T) 
i 7t-7l, (l=(l(T) 

7i-h, (i=(i(T) 



FIG. 5: Relaxation time of shear viscous flow as a function of temperature 


temperature range shown. When /i = /i(T) is used, the density n is larger at higher tem¬ 
peratures where fi is smaller. Consequently, this compensating effect is enhanced resulting 
in an almost insignihcant variation of with temperature. 

Note that the relaxation times for dissipative flows which have been plotted in hgs. [3l H] 
and [5] on the same scale come out to be of similar magnitude. The hrst order coefficients 
however are quite largely separated. The bulk viscosity ( is generally much smaller than rj 
as seen in e.g. [l^, 15 1. Consequently, bulk viscosity and thermal conductivity are usually 
ignored in the set of hydrodynamic equations. 


V. SUMMARY 

In this work we have evaluated the relaxation times of dissipative fluxes in the kinetic 
theory approach using Grad’s 14-moment method for the case of a pion gas. Our aim has 
been to estimate the in-medium effects on the temperature dependence of the relaxation 
times. Using an effective Lagrangian the tttt scattering amplitudes involving p and a meson 
exchange have been evaluated in which one-loop self-energy corrections were incorporated to 
obtain the in-medium propagators. The consequent decrease in the effective cross-section is 
found to have an appreciable effect on the temperature dependence of the relaxation times 
of the irreversible flows. Since these go as inputs in the second order viscous hydrodynamic 
equations it is expected that the space-time evolution of heavy ion collisions will be affected 
signihcantly. So a realistic evaluation of these quantities is essential to obtain the proper 
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temperature profile and consequently the cooling laws of the evolving system. In addition it 
is found that the relaxation times of the bulk viscous flow and the heat flow to be of similar 
magnitude to that of the shear viscous flow which suggests that they should all be taken 
into consideration in dissipative hydrodynamic simulations. 


Appendix A 


The equations of motion of the thermodynamic variables follow from the conservation of 
particle number and energy-momentum along with the contraction of the last quantity with 
the hydrodynamic four-velocity and the projection operator The evolution equations 
of particle number, four-velocity and energy density along with those of temperature and 
chemical potential for a dissipative fluid are respectively given by 


Dn = —nd ■ u. 
hnDu>^ = V^P - AyM’"'' - 

nDe = -Pd ■ u - 
T-^DT = (1 - y)[9 • n + 

TD{^^} = {(i' - l)h - i"T]\d . V. + - n“-V.u„}] 

--{V.y - n»“'v„u„}. 

n ^ 

where, 

^ (S»/S‘)2 - + 32-iS»S‘/(S‘)2 + 2z->S«/Si - z-^- 

2-2 

^ - (SI/Sir- + 3*-'S»S‘/(S>)2 + 2z-iS»/S‘ - 2-2 

_ Sj7s^ + 52-"s|/s^ - sisy(sif 

^ (S?/S‘)2 - (S»/S72 + 3z-iS»SP(S‘)2 + 2z-^SI/S^ - z-^ 


(55) 

(56) 

(57) 

(58) 


(59) 

(60) 
(61) 
(62) 


s|S27(S2‘)^ 

1 - nsssi - Ss‘S2”}/(Si)" 
-1 


(63) 

(64) 


In the first order theory all the flows and dissipative quantities are ignored in the above 


equations of motion since these are first order in the gradients, while the thermodynamic 
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variables are of zeroth order 


20| . So in Chapman-Enskog method the equations of mo¬ 


tion fl551l59p have been used as thermodynamic identities excluding the dissipative terms 
like and fl^^" as given in |l^. In case of the second order theory the space of thermo¬ 
dynamic quantities is expanded to include the dissipative quantities which are treated as 
thermodynamic variables in their own right. However, close to equilibrium the gradients of 
macroscopic variables may be treated as small quantities of the order of 0. Also, terms like 
and which are products of the fluxes and gradients are of higher order in 

(f). Following 2^ we neglect such higher order terms in equations 


Appendix B 

The zeroth order distribution function contains a few arbitrary parameters which 
are identihed with the temperature, chemical potential and hydrodynamic 4-velocity of the 
system by asserting that the number density, energy density and the hydrodynamic velocity 
are completely determined by the equilibrium distribution function. We thus write the 
number and energy densities as 

n = gj Tpp^uJ = g j (65) 

en = g J rp(p%)V = 9 J (66) 

where g = 3 is the isospin degeneracy of pions. Furthermore, from Fckart’s dehnition of the 
hydrodynamic velocity we get. 


= j = 0 . 


(67) 


From the above three equations we obtain equations ffT^ . fl^U]) and fl^ which on substi¬ 
tution in the expression for 0 results in. 


1 r,rA(x,r)/(°)(l + /(°)) = 0, 

(68) 

1 r,r2A(x,r)/(°)(l + /^°^) = 0, 

(69) 

1 r,(n05.(x,r)(n^)/(°)(i-t-/(°)) = 0. 

(70) 


Expanding the coefficients A and according to equations ([9]) and (fTOj) . we hnally 
obtain ([22]), fl23l) and (1241) . 
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Appendix C 


Here we elaborate on the coefficients of bnlk viscosity, thermal conductivity and shear 
viscosity which are given by 


c 

A 

V 



022 


T P! 

3m^ bii 
10 Coo 


(71) 


The numerators have been dehned earlier. The quantities 022 , &11 and coo are dehned in 
terms of the collision bracket 


[F,G] = ^J dTpdTkdTp, dr,,/(°)(p)/(°)(fc){l + /(°)(p')}{l + /^°Hfc')}5(^)5(G) W (72) 
with 


In particular, 


d(F) = F(p) + F(k) - F(p') - F(k') 
S(G) = G(p) + G(k) - G(p') - G(k') 


“22= IAt% 

bn= -A^„[T{n'‘).r(n‘')], 

m 

coo= ^[(n^n,),(n^n")]. 

Tflj 

They can be expressed in terms of the integral X^iz) as 

O22 = z‘^X^{z) 
bn =-z[X2{z) + Xs{z)] 


and 


where 


liml 


Coo — 2 [JVi(c) + -^2(0) + -A3 (c)] 


(73) 

(74) 
(76) 

(76) 

(77) 

(78) 

(79) 


XJz) = 


Sz* 


|S](c)P 


,{-W) 


/•TT pZTT 

difj cosh^ -0 sinh'^ -0 / d© sin 0 - —( 0 , 0 ) / dcj) 
Jo zaiZ Jo 


dx sinh^“ X I d6 sin 6 


^2z cosh ijj cosh x 


(e^-l)(e^-l)(eG-l)(e^-l) 


M,(d, 0 ) (80) 
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with 


E = ; 2 (coshcosh X — sinh'^sinhx cos 6^) — /i/T 
F = ; 2 (coshcosh X — sinh'^sinhx cos 0') — /i/T 
G = E + 2z sinh ip sinh x cos 9 
H = F + 2z sinh ip sinh x cos 9' . 


The functions stand for 

Mi(6*, 0) = 1 — cos^ 0 , 

M 2 (0, 0) = cos^ 9 + cos^ 9' — 2 cos 9 cos 9' cos 0 , 
M3{9,Q) = [cos^ 0 — cos^ 6*']^ 


and 


cos 9' = cos 9 cos 0 — sin 6* sin 0 cos (p . 


(81) 


(82) 


(83) 


Appendix D 


In this appendix we provide explicit forms of some integrals which have appeared in the 
expressions for the transport coefficients in the text. Recall that an = f 
from which we obtain 


Oi = 


02 — 

03 = 

04 = 


05 — 


n ^ 


1 }, 






n 


C2 

-2 




n Ql Q3 cO c2 

+ {^ + 

1 O2 O2 O2 O2 


(84) 
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and from = f drpf(^^(l + /(°))r"'(n^)(n'^) we have 


f. ^ 

0o = - 

h =_ 

T^Sl’ 

62 = - ^{ 5 z ^ + z ‘^}, 

, ^ron >^1 r 2>^l 3 >^ 3 S 

h — — 77;{302:-^ + 5z -^ + z -^}. 

i O2 ^2 *^2 

Again, = -^ / drp/W(l + /®)(2r’" gives 

1 cO QO q cl Cl 

«2 = ;.^[-(^ - ;.-^) + (^ + -^){(1 - 7 ")^ + rz -^)} 

O O2 *^2 ^ ^2 *^2 

-4-7'){| + 15.-^|| + 2 .-'}|. 

O O2 *->2 

/^n = / drp/(°)(l + ■ M - h)Ap^p^^p'' gives 

ft = 3*-ni + 5.-‘|-4n. 

*^2 

and from 7 ^ = / dTpf^^\l + f^^^)r'^{ppPu){p^p'') we get 


c2 

■'»- -“4 

C3 c2 

71= -10{^ + ;^f}. 

*^2 *^2 


(85) 


( 86 ) 


(87) 


( 88 ) 
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